Phase velocity and phase diffusion in periodically driven discrete state systems 
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We develop a theory to calculate the effective phase diffusion coefficient and the mean phase 
velocity in periodically driven stochastic models with two discrete states. This theory is applied 
to a dichotomically driven Markovian two state system Explicit expressions for the mean phase 
velocity, the effective phase diffusion coefficient and the Peclet number are analytically calculated. 
The latter shows as a measure of phase-coherence forced synchronization of the stochastic system 
with respect to the periodic driving. In a second step the theory is applied to a non Markovian 
two state model modeling excitable systems. The results prove again stochastic synchronization to 
the periodic driving and are in good agreement with simulations of a stochastic FitzHugh-Nagumo 
system. 
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Stochastic resonance as a phenomenon of noise en- 
hanced order in periodically driven stochastic systems 
has attracted considerable interest until today [H El 0- 0| . 
A common approach to quantify this effect are spectral 
based measures like the spectral power amplification and 
the signal to noise ratio. On the other hand stochastic 
resonance can also be understood as a synchronization 
process between the input and the response of the sys- 
tem @, l|- This interpretation achieves importance es- 
pecially if dealing with larger amplitudes of the driving 
signal. Then analytical descriptions have to go beyond 
linear response theory. 

In general two principal approaches were introduced 
in the past to describe the synchronization of a stochas- 
tic system by an external driving. The first one bases 
on the consideration of escape time densities to leave 
certain states of the dynamical system 0. A periodic 
driving modulates these densities and they exhibit max- 
ima at times which corresponds to time scales of the ex- 
ternal drive. "Bona fide" resonances were investigated 
analytically, numerically simulated and and experimen- 
tally verified, especially for symmetric bistable situations 

The second approach goes back to Stratonovich who 
looked at synchronization of nonlinear oscillators by pe- 
riodic driving in the presence of noise ^lj • F° r this pur- 
pose one adopts a phase to the nonlinear oscillators and 
defines statistical properties of the stochastically behav- 
ing phase. If the mean phase velocity agrees with the 
frequency of the driving and at the same time the phase 
diffusion coefficient is small then there exist in average a 
fixed phase relations between the driving and the output 
of the system. 

This picture was recently transfered to models of 
stochastic resonance which are nonlinear but non os- 
cillating. It was possible to prescribe a phase to over- 
damped bistable as well as to excitable systems which 



monotonously increases in time EHHHl- Its mean ve- 
locities and effective phase diffusion constant were used 
to quantify synchronization between the output and the 
driving input. Likewise as in stochastic resonance syn- 
chronization appears at an optimal choice of the noise 
intensity since the level of noise determines the charac- 
teristic times of the stochastic system. 

As result one finds plateaus of the mean frequencies 
of the output at values which corres pon d to the driving 
frequency or multiples of it 0, EH [01 E3, These 
plateaus are accompanied with low phase diffusion coef- 
ficients indicating a synchronization in average. As mea- 
sure of synchronization one uses the duration of lock- 
ing epoches or a Peclet-number which is the ratio be- 
tween the phase velocity and phase diffusion coefficient 

B El Lull . 

For bistable stochastic systems a discrete state model- 
ing has been proven very successful in the past (2^. l2^|. 
It is based on a separation of time scales between the fast 
relaxation into the metastable states and the transition 
between these states, which happens on a slower time 
scale and build up a Markovian discrete dynamics [2J|. 

Also models of excitable behavior [2E 15(1 can be 
mapped on two or three state dynamics jU EE HH- 
These discrete state models still set up a renewal pro- 
cess However in difference to bistable systems they 
include non-exponentially distributed waiting time den- 
sities and are thus non Markovian. 

These discrete state systems will be endowed with a 
discrete phase which is introduced in Section [H] As will 
be shown in our paper both the Markovian and the non 
Markovian model exhibit phase synchronization with re- 
spect to the periodic driving for optimal noise levels. We 
will quantify this effect by the mean phase velocity, phase 
diffusion coefficient and the Peclet number. An unique 
approach to calculate these quantities in driven renewal 
models with two states will be presented in Section HTT1 
This a ppr oach is based on an envelope description of the 
phase HE EH- 
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Section llVl anplies the theory to bistable systems where 
Markovian rules were assumed for the transition between 
the discrete states. First results of this system with 
dichotomic periodic inputs were derived earlier in [lflj . 
These results were recently improved in [32I l8^ | which 
agrees with our findings in case of Markovian dynamics. 

Section is devoted to a non Markovian two state 
system which models excitable behavior. Integral equa- 
tions for the phase velocity and phase diffusion coefficient 
have to be numerically solved. Results of these compu- 
tations show good quantitative agreement with numeric 
simulations of a stochastic periodically driven FitzHugh- 
Nagumo system. 



II. TWO STATE MODELS AND PHASE 

Consider a periodically driven stochastic two state sys- 
tem described by the probabilities p(t) = \pi(t) , to 
be in state 1 or 2 respectively at time t. 

In generally these dynamics can be expressed in terms 

by 



of the flux operators Jl~^ 3 



Pi 

I>2 



2r P(-)] 



+2 [p(-)i 
+i [p(-)i 



(i) 

(2) 



The linear flux operators, which express the probability 
flux from state i to state j in terms of the occupation 
probabilities depend explicitly on time t in a periodic 
way due to the periodical driving with period T — 2-k/SI, 



Jt 



(3) 



In the Markovian case these operators are local in time, 
i.e. multiplication operators, 

The well known two state model for bistable systems 
which will be considered in more detail in section Hvl is 
of this type. In the non Markovian case the action of 
the flux operators Jl~^ 2 on the probabilities p\ and P2 is 
non-local in time, i.e. the are integral operators. 

One example of this type is the discrete state model for 
excitable systems [2Sl. Efl| , whose flux operators are given 

by 

J^ 2 [p{-)]{t) = / d™(t-r) 7 (r)pa(r) 

Jt 

jr^vi-m = 7(*)j*(t) 

Note that in this case the flux operator depend explicitly 
on the initial time to which breaks its periodicity eq. {3J. 
However in the asymptotic case to — » —00 this periodicity 
is restored. This model will be considered in section Fvl 

Next we endow this system with a phase <fi(t). Our 
goal is to evaluate the mean phase velocity 



uj := lim 

t — >oo 



(4) 



as well as the effective phase diffusion constant 



D, 



eff 



:= lim 

t — >oo 



2t 



(5) 



These quantities are independent of the exact definition 
of phase, as long as the phase increases by 2ir within a 
one cycle 1 — > 2 — * 1 of the system. For the sake of 
notational and computational convenience we consider a 
phase, which increases by 2ir each time the system enters 
state 1. Then the probabilities Pk — [pi,fc 5 P2,fc] to be in 
state 1 or 2 respectively and to have the phase 2nk are 
governed by 



Pi.k = J^ l [Pk-i]-J^ 2 [Pk] 

P2,k = Jt^iPk}- Jt^\Pk] 



(6) 
(7) 

These equations are similar to eqs. Q and 0, however 
the probability influx into state 1 for a given phase 2itk 
comes now from states with the phase 27r(fc — 1). 

The mean phase as well as the mean square phase are 
given in terms of the probabilities pfe by 

00 

k— — oo 

00 

W*)) = E ^ 2 k 2 (piAt)+P2Mt)) 

k— — oo 

The instantaneous mean phase velocity ui(t) and instan- 
taneous mean phase diffusion D e g(t) are then defined as 



w(t) 
D eS (t) 



dt 
1 d r 
2dt 



W)) 



(8) 
(9) 



Asymptotically, i.e. for the initial time to - * —00, the 
phase (j> = 2irk will undergo a diffusional motion with 
periodically varying effective phase velocity tj(t) and ef- 
fective diffusion coefficient Z? e fF(^)- In this asymptotic 
regime the mean phase velocity eq. I@J and effective 
phase diffusion constant eq. © can be expressed as the 
time average over one period of the external driving of 
the time dependent phase velocity and diffusion constant, 

Q = i / dtuj{t) and D eS = i / dtD eS (t) (10) 



T 



T 



Although the phase velocity and effective phase diffusion 
constant eqs. JBJ and JjjJ have a periodic asymptotic 
behavior, the equations (jSJ and Q which govern the 
probabilities p\± and p2,k obviously have no asymptotic 
solutions. 



III. GENERAL THEORY 

Our aim is to relate the asymptotic phase velocity and 
effective phase diffusion constant eqs. ifTTHl to the mi- 
croscopic dynamics eqs. Q and 0. To this end we 
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introduce a continuous phase distribution V(4>, t) as the 
envelope of the discrete phase distribution pi.k and p2,k 
[20L 0] by defining its values at integer multiples of 2n 



V{(f) = 2nk,t) :=Pi,fc(t)+P2, fc (t). 



(11) 



The diffusional motion of the phase <j> requires its distri- 
bution V(4>) to obey the Fokker-Planck equation 



0_ 



n*, *) = + D eS (t)^)p(<p, t). (12) 



To establish the relation between oj(t) and D e g(t) and 
the microscopic dynamics eqs.l© and (0 we expand pi : k 
and p2,k according to 



<~*~> fin 
n=0 V 



-2wk 



1,2 (13) 



This expansion describes how the probability to be in 
state 1 or 2 for a given phase k at time t, Pi,k(t) an d 
P2,k{t) respectively, is related to the total probability to 
have a phase 2irk, T(2itk,t) and its gradients. 

The total probability Pi,k{t) + P2.k(t) to have a phase 
2irk neglecting the internal state 1 or 2 is related to the 
continuous phase distribution by the defining eq. lilt , 
which in turn implies 



(*) = 



(0; 



(14) 



for n > 1. 



Inserting the Ansatz eq. 113t into the master eq. (j^J and 
Q , using the Fokker-Planck equation 11211 for the phase 
and considering the coefficients of the different deriva- 
tives d n /d4> n V{4>, t) eventually leads to (cf. appendixEJ) 



q(°) = M t [q (0) ] (15) 
q« = Mt[q( 1 )+c| 1) q^]-27rJ I t n [q^]+w(-)q (0) (16) 
q< 2 ) = M 4 [q (2) + C^qW + cf\^\ - 27rJ»[qW + (c^ - n)^] + c(-)q (1) - D eff (-)q^ (17) 



where we have introduced the master operator 
The operator 

Jf H = ( J \ 1[ ' ] ) 
accounts for the influx into state 1 and we introduced 
cfV) = / drua(r) 

i-t r-t 

cf\t') = - J cLtD {t) + wo(t) J t dr(t - t)wo(t). 

q(°' in eq. l(T5)) shows the same dynamics as p in the 
two state system without phase eqs. and |2"| l. which 
one would also expect as this term corresponds to an 
equipartition of phases t) = const in the expansion 
eq. 11 31 . The higher order terms q(") are corrections 
which emerge due to the fact that we are considering a 
non equipartition of phases resulting in drift and diffu- 
sion. 

Interestingly if the action of the flux operators on the 
probabilities is local in time, i.e. in the Markovian case 
the terms containing the Cj are zero, as Cj (t) = 0, and 
therefore the dynamics of the q^ considerably simplifies. 



By summing up both components of the vectorial eqs. 
itTfi)) and itTTjl , using the normalization condition eq. ifHjl 
and the fact that (M t )i + (M f )2 = we arrive at 

W (t) - 2 7 rj^ 1 [q (0) ]M (is) 

D eS (t) = 2 7 rJ 4 2 - 1 [-q (1) + (^-ci 1) )q (0) ]W 

= nu(t) + 2<7 t 2 ^ 1 [-q (1) - c'^q^Kt) (19) 

The asymptotic mean phase velocity Q and the asymp- 
totic effective phase diffusion constant D e g can then be 
determined from the asymptotic (cyclo stationary) solu- 
tions of eqs. I|15|) and 11 (ill . Therefore, the calculation 
of the asymptotic effective diffusion constant is reduced 
to the solution of a cyclo stationary problem, which in 
general is simpler that solving the whole non stationary 
equations 10 and {JJ with some initial conditions and 
then taking the asymptotic limit in eq. ©. 

In the following the mean phase velocity and effective 
phase diffusion constant will be considered for two differ- 
ent models, namely a Markovian model [23], which ap- 
proximates bistable systems and a non-Markovian model 
|29l |. which serves as an approximate description for ex- 
citable systems. For the dichotomically driven Marko- 
vian case the mean phase velocity and effective phase 
diffusion constant can be explicitly calculated, while for 
the non-Markovian case solutions can only be obtained 
numerically. 
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IV. A MARKOVIAN TWO STATE MODEL 

We consider now a Markovian two state system with 
periodically modulated rates 72 (0 and 71 (t). Its flux 
operator J 1 ^ 2 and J 2 ^ 1 are given by 

#-*[!>](«) = 7i(*)Pi(*) and ^[pK*) - Ts(t)Pa(t) 

In this Markovian case, the equations, which govern the 
evolution of qW (£) greatly simplify due to the fact that 

cjj* } (i) = 0. Eqs. © and JEJ reduce to 

w(t) = 2^ 72 (t)g( 0) (t) 
£>eff(*) = 27r 2 72(*)92 0) W- 2^72(4)9^ (t). 



The equations for an d #2 (*) are gi ven by 



(0), 



(20) 



= 7iW9i 1) W-72W^ 1) W+^W4° ) (*) (21) 

Eqs. (f2T)j> and 1I2TI) can be readily solved by the method 
of variation of constants, using qf\t) = 1 — q^ 2 °\t) an d 

ffi^C*) = — 9*2 (*) ( c f- e 1- CU)- The asymptotic periodic 
solutions eventually read 



Jo dT H (t - T ) exp(-s(r, t)) 
l-exp(-s(T,*)) 



(22) 



(l),., Jo dTUJ (t - T)qf ] {t - r) exp(-s(r, t)) 
q2 W ~ l-exp(- S (T,i)) (23j 

where s(r,t) := f*_ T dr' (71 (r') + 72 (t'))- Note that 
s(T,t) does no longer depend on t. 

For a dichotomic symmetric driving with period T = 

2tt/w, 



, , n if ie[nT,(n+±)T) 

1 r 2 if te [(n + i)T,(n + l)T) 



and vice versa for 72 (t) eqs. 122t and 12 -ill can be readily 
evaluated leading after some cumbersome algebra to the 
mean phase velocity and effective phase diffusion con- 
stant 



£1 tanh R 



and 



D e e = 7rw [i +a(i +cosh 2 R)] 

,1 



(24) 



(25) 



7rafitanhi?[ - 1 + a(- cosh 2 i? + 1)] 

where we have introduced the mean phase velocity with- 
out driving wo := 2w/(^-+^), a quantifier for the driving 

strength a = and some ratio between inner time 

scale and driving frequency R = ^-^f^-- 



Without signal, i.e. a = eq. 112511 reduces to 
-D e ff = Tr^o, which agrees with the result in 
(2^) 2 /2«i 2 )-(i) 2 )/(t) 3 . 

Next we consider the small and large noise limits of the 
phase velocity u> and phase diffusion constant D eff for the 
case of Arrhenius rates ri/2 = exp(— AU fj zA )- In this 



case a — tanh 2 (-^). 

If for a fixed driving frequency the noise level is suffi- 
ciently small such that R -C 1 eqs. l(2"4*jl and (H} reduce 
to 

7T 7T 

uj w uj + aQR = -(n + r 2 ) w -r 2 

13 3 
D eS w 7tw (- + -a) + 7rafti?(-l + -a) 

7T 2 . . 7T 2 

= T (ri+r 2 )« T r 2 

where in the last step we used the fact that r 2 dominates 
7"i for small noise levels. Therefore, at the level of phase 
velocity and phase diffusion, the process behaves like a 
process without driving whose rates are both equal to I j. 

On the other hand if the noise level is large and the 
driving frequency is small compared to ag such that R ^> 
1 we get 



r\r 2 
n + r 2 



n 



(ri - r 2 ) 2 



-ui (l + a)+ 7raO(-l + a) 



2tt 



; rir 2 (r 2 + rg) 
(ri + r 2 ) 3 



4tt^ 



(ri + r 2 ) 2 

(n - r 2 ) 2 
'(n + r 2 ) 4 



The first terms in these expressions correspond to a pro- 
cess without driving with one rate equal to n and the 
other equal to r 2 , while the second terms are corrections 
which vanish for vanishing driving frequency. 

Between these regions we have a competing behavior. 
If for a fixed driving amplitude A, the noise strength 
D is sufficiently small, such that awl and ujo w 0, 
and simultaneously, for a fixed driving frequency fi, D is 
sufficiently large such that R 3> 1, i.e. tanhi? w 1 we 
have 







i.e. frequency and phase locking occur. 

Having calculated the effective diffusion coefficient and 
the mean phase velocity we can evaluate the Peclet num- 
ber 



Pe := 



2ttuj 



D 



eff 



(26) 



which is a measure of the phase coherence. 

In FigQ]the theoretical results eqs. lj24l) . ll25)l and J2£ 
are compared to simulations of the driven two state sys- 
tem. To compute these results we have modified an al- 
gorithm presented in j^^l taking into account that the 
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transition rates are piecewise constant in time due to 
the dichotomic driving. Let us assume we start at time 
t in state 1 and the input defines the rate to have 
the value r\. Then we draw a random number t ac- 
cording to the corresponding waiting time distribution 
W ri (r) = r\ exp(— t\t). If t+r is smaller then the time t s 
of the next switching of the input we set the running time 
to t + t and perform the transition to the second state 
of the system. This state 2 will be left with rate r 2 and 
we proceed accordingly. Contrary if during the interval 
[t 7 t+r] a switching of the input occurs we set the running 
time equal to the switching time t s but remain in state 1. 
After switching of the input the rate for leaving state 1 
is now r 2 and we proceed by drawing a new waiting time 
according to the new density to r2 (r) = r 2 cxp(— r 2 r). 




0.005 0.01 0.015 0.02 0.025 0.03 0.035 0.04 
D 



Figure 1: Mean phase velocity u)q (top), effective phase diffu- 
sion constant Do (middle) and Peclet number Pe (bottom) of 
the Markovian model for different values of the driving am- 
plitude. 

Symbols are simulation data of the two state system, lines 
according to eq. 12411 . 125t and 12611 . respectively. Other pa- 
rameters: ro = 1 and AU = 0.25, fl = 0.0017T. The deviation 
between theory and simulations in the Peclet number for low 
noise intensities is due to limited simulation time. 



of noise strength, indicating stochastic resonance. For a 
strong driving, it varies over several orders of magnitude 
with varying noise strength D. Interestingly the Peclet 
number shows also a non monotonic behavior as a func- 
tion of driving frequency for a fixed noise level, i.e. using 
this number as a measure of the quality of the response 
to the external signal we discover a "bona fide " resonance 
(Fig.©. 




10' 1 1 ' 

0.01 0.02 0.03 0.04 0.05 0.06 

£2 



Figure 2: Peclet number Pe of the Markovian model as a 
function of driving frequency f2 for different noise values. A = 
0.2, other parameters as in Fig. Q The inset shows the driving 
frequency at the maximum Peclet number as a function of 
noise strength (solid line) compared to the intrinsic frequency 
r without driving (A = 0) (dashed line). 



V. EXCITABLE SYSTEMS 

In this section we consider the phase velocity and dif- 
fusion of a non Markovian model [2^|. This two state 
model mimics the dynamics of an excitable systems by 
dividing it into an excitation step and the evolution along 
the excitation loop. Its dynamics is given by 

pi(i) = 7(*M*)- / drw(t - t)j(t) P2 (t) (27) 

Jto 

p 2 (t) = - 7 (*) 1 fc(t)+ f d™(t-r) 7 (r)p 2 (r). (28) 

Jto 

with initial conditions 

pi(to) = and P2 (t Q ) = l. (29) 

State 2 represents the rest state, in which we start at 
initial time to >From there the system is excited due 
to noise and the external periodic subthreshold signal, 
leading a rate process with rate 7(i), which depends pe- 
riodically on time. This Markovian excitation step is 
described by 



The Peclet number shows a maximum as a function 



^ x tp](*) = 7(«*)- (30) 
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State 1 accounts for the motion on the excitation loop on 
which the systems spends a time distributed according to 
the waiting time distribution w(t), which is assumed not 
to depend on the weak external driving. The flux from 
state 1 back to state 2 is then expressed in terms of the 
flux from state 2 to state 1 at prior times r between to to t 
i(t)p2(t), which renders the description non Markovian, 
leading to the flux operator 

J t ^ 2 [p}(t) = [ drw(t-Th(T)p 2 (r). (31) 

Note that this operator depends explicitly on the initial 
time to. 

To calculate the asymptotic periodic solution it will be 
useful to first formally integrate eqs. ffijl . jSEt taking 
into account the initial conditions H29H and then taking 
the initial time to to — oo. The resulting equations are 

roc 

Pl (t) = / drz(T)j(t - r)p 2 {t - t) (32) 
Jo 

/>oo 

P2 (t) = 1- dTZ(T)f(t-T) P2 (t-T). (33) 

JO 

where z(t) = 1 — J Q T dr'w(T') is the probability to spent 
a time longer than r on the excitation loop. By differen- 
tiating these equations with respect to t one recovers the 
original eqs. iJSZl and <GBJ in the limit t — > -oo [35] . 

If we take into account the phase eq. l(33jl has to be 
replaced by 

/>oo 

P(i,k)(t) = / rfT^(T)7( i -' r )P(2,fc-i)(i-T)- (34) 
Jo 

We also have to take care of the flux operator j}^ 2 which 
in the asymptotic case is given by (cf. eq. {BU) 

/•OO 

^ 2 [p](i) - / drw{T)j(t-r) P 2{t-r). (35) 
Jo 

In the following we assume a fixed waiting time T on 
the excitation loop, i.e. w(t) = S(T — t) and z(t) = 
8(T — t). Such an assumption is justified in the low noise 
limit for e.g. FitzHugh-Nagumo models (cf. Fig. UjJ. In 
this case eq. i|35ll simplifies to 

J t ^ 2 [p](t) = y(t-T)p 2 (t-T). (36) 

Then, according to eqs. iflSll and iflfljl . the time de- 
pendent phase velocity ujq (t) and effective phase diffusion 
constant Doit) are given by 

w(t) - 2^{t) q { 2 a \t) (37) 
D eS (t) = ~2ir 1 (t)qi 1) (t) + 2ir 2 1 (t)qi a) (t), 

which are the same expressions as in the Markovian case, 
as the flux operator J 2 ^ 1 is the same. However the 
equations governing the qW are different. Following the 



same procedure we used to treat eqs. 10 and eq. l.'i-lll 
together with normalization the condition itTHl leads to 

l-$\t) = [ T dr 7 (t-T)q^(t-T) (38) 
Jo 

-<£\t) = f d T1 {t~T) q { i\t-T) (39) 
Jo 

+ d T1 {t - T)q ( 2 ] (t - t){ j dT'uj{t-T') -2n). 
Jo ~ Jo 

The periodic solutions of eqs. il38|l . l(39"jl can be numeri- 
cally obtained in Fourier space using a linear solver like 
LAPACK. 

To investigate the role of noise on the synchroniza- 
tion in excitable system we choose an Arrhenius type 
excitation rate for the transition from the rest state 2 
onto the excitation loop 1. We further assume that the 
external driving acts as a modulation of the potential 
barrier. Again we consider a dichotomic periodic driv- 
ing, i.e. the excitation rate ^(t) periodically switches 
between the two values r\ = ro exp(— (AU — A)/D) and 
r 2 =r exp(-(A[/ + J 4)/D). 

The resulting phase velocity, effective phase diffusion 
and Peclet number as a function of noise strength D are 
shown in Fig. As in the case of bistable systems we 
observe frequency and phase locking, however there ex- 
ist preferred driving frequencies for which high synchro- 
nization is achieved and other frequencies which show no 
synchronization at all. 

The Peclet number shows a local maximum at a finite 
noise strength. Contrary to the bistable situation how- 
ever, the phase diffusion constant decreases again and the 
Peclet number therefore increases for large noise levels. 
This behavior is originated in the fixed time T on the 
excitation loop. Taking into account the high rate and 
therefore small waiting time and variance of the excita- 
tion step for high noise levels this leads to a low variance 
of the spiking, which implies a low diffusion of the phase. 
We mention that this low phase diffusion does not im- 
ply synchronization since the frequencies are not locked. 
Also we note that in real excitable systems the behavior 
differs. For higher noise levels the time spent on the ex- 
citation loop will have a variance in these systems which 
yields an increasing phase diffusion with growing noise. 

As seen in Fig0the synchronization behavior strongly 
depends on the driving frequency. To further analyze 
this effect we have plotted in FigQ]the mean phase veloc- 
ity, phase diffusion coefficient and the Peclet number as 
function of the driving frequency. They show a complex 
sequence of different locki ng regio ns between the driving 
and the system's response [36Ll37l |. represented by shaded 
regions. In these locking regions the effective phase dif- 
fusion is small (see FigEJ ■ We mention that the maximal 
frequency of the excitable system is c^max = 2n/T where 
T is the time on the excitation loop. There can not be 
1 : 1 synchronization for £1 > ciJ max . 

Let us for a moment assume the extreme case where 
one excitation rate n is infinity and the other r 2 is zero. 
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Figure 3: Mean phase velocity Q (top, inset), effective phase 
diffusion constant D e g (top) and Peclet number Pe (bottom) 
of the non Markovian model for different values of the driving 
frequency SI. 

Symbols are simulation data of the two state system, lines 
according to numerical evaluation of the theory. Other pa- 



rameters: T - 
A = 5.0- 10" 5 . 



2800, ro = 0.0044, AU = 5.6 • 10" a and 



Figure 4: Mean phase velocity Q (top), effective phase diffu- 
sion constant D e g (middle) and Peclet number Pe (bottom) of 
the non Markovian model as a function of driving frequency 
Q, for D = 0.00001. The shaded regions are a guide for the 
eye and represent regions of frequency synchronization. In 
these regions we also find a small effective phase diffusion and 
therefore a high Peclet number. Other parameters: T = 2800, 
r = 0.0044, A 17 = 5.6 • 10 -6 and A = 5.0 • 10 -6 . 



Then the system remains in the rest state as long as 
the input causes the vanishing excitation rate. After the 
input changes the system immediately starts with the 
excitation loop where it stays the time T. For a 1 : 1 
locking this time T must be larger than half the period 
but smaller the full period 27r/f2 of the driving. Oth- 
erwise, if the duration of the excitation loop would be 
smaller than half the period the system returns to the 
rest state where it immediately starts a new excitation. 
In consequence 1 : n locking where the output frequency 
is n times higher than the input frequency occurs if the 
period of the driving is between (n — 1/2)T and nT . 

The opposite case where a fast input locks a slow out- 
put occur if multiple periods of the input fit into the 
excitation time. During the excitation the system does 
not respond to the changes of the input. If the input has 
the phase with long waiting time after the system has 



completed the excitation loop, it has to wait until the 
input changes to the phase with the small waiting time, 
leading to a n : 1 synchronization where n is the number 
of signal periods which fit into the excitation time T. 

However if the system finds the high excitation rate 
after excursion it immediately starts a new excitation 
loop and repeats these until it will find the phase with 
long waiting times. This yields a n : m frequency locking 
with n > to. Note that there are no n : to locking modes 
with n < to except the 1 : to modes described above. 

Realistic noise dependent time scales will weaken the 
extreme behavior of the situation considered above. 
There are two competing effects namely increasing the 
noise increases r\ as well as r-i while decreasing the noise 
increases the ratio between r\ and r 2 and therefore the 
effect of the driving. Hence, we find synchronization in a 
finite window of noise intensities where the two activation 
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Figure 5: Effective phase diffusion constant D e g of the non 
Markovian model, as a function of driving frequency f2 and 
noise level D. The black lines show regions of frequency lock- 
ing (1 - e)fi < nui < (1 + e)fi, e = 0.01 with (from left to 
right) n — 3,2, 1, |. These regions of frequency locking coin- 
cide with low phase diffusion. Other parameters as in Figs. El 
andU 



times enclose the time T on the excitation loop, 

1 1 
— <T< — . 

n r 2 



(40) 



We point out that this latter time plays the essential role 
within the synchronization process, i.e. this time scale 
and the period of the external drive have to be tuned 
appropriately to get phase synchronization. Noise as well 
as the amplitude of driving define the two excitation rates 
and have to be chosen such that eq. lUTHl is optimally 
fulfilled, i.e. that the input acts as much as possible as a 
on-off switch on the excitation process. A deviation from 
this extremal behavior leads to a narrowing of the driving 
frequency windows amenable to frequency locking and a 
shift of these windows to lower frequencies. 

Finally we compare the theory to a dynamical system 
with excitable dynamics, namely the FHN model 
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Figure 6: Inter spike interval distribution of the FHN system 
eqs. I4ip with constant signal s(t) = 0.04 for a low noise level 
D = 10~ J and strong time scale separation e = 0.001. Other 
parameters see text. 




Figure 7: Comparison between theory and FHN system for 
mean phase velocity uj (top) and effective phase diffusion con- 
stant Z) e ff (bottom) as a function of driving frequency. Other 
parameters see text. 



X = x - x A - y + V2D£(t) 
V = e(x + a - aiy - s(t)) 



(41) 



This system is driven by a dichotomic periodic signal s(t) 
with values ±A where A = 0.015. Setting ag = 0.405 and 
CH = 0.5 the system is in the excitable regime for both 
values of the signal, i.e. the signal is a sub-threshold 
signal. We further consider a strong time scale separation 
e = 0.001 as well as a small noise level D = 10 -5 . The 
phase of the system is defined to increase by 2ir each time 
a spike is generated. From simulations of the inter spike 
interval distribution (see Fig. EjJ for constant signal ±A 



we find the corresponding parameters of the two state 
model to be T w 2620, n « 0.0087 and r 2 « 8.310~ 8 . 

The results for the phase velocity Q and effective phase 
diffusion constant D e g for the FHN system (numerical 
simulation of eqs. ijllll ) and the theory eqs. ll3Tjl and 
l|38ll are shown in Fig. They show a good qualita- 
tive agreement over a large range of driving frequencies. 
The deviation for larger driving frequencies is due to the 
fact that, in contrast to the assumptions of our two state 
model, the time T spent on the excitation loop depends 
if however only weakly on the driving. 
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VI. CONCLUSIONS 

We have derived a general theory to calculate the 
asymptotic effective phase velocity and phase diffusion 
constant in periodically driven two state systems. This 
theory was applied to two different two state models, one 
with Markovian dynamics representing bistable systems 
and the other with non Markovian dynamics, modeling 
excitability. 

In the Markovian case analytical results have been cal- 
culated for dichotomic driving with arbitrary driving am- 
plitudes. We found phase synchronization for optimal 
noise intensities if Arrhenius type rates for the transitions 
between the states are assumed. The mean frequency of 
the system is locked to the frequency of the external stim- 
ulus and the effective phase diffusion coefficient becomes 
vanishingly small. The Peclet number however shows a 
maximum not only as a function of noise strength but 
also as a function of driving frequency, i.e. a "bona fide" 
resonance. Frequency locking occurs as long as the driv- 
ing frequency is smaller than the maximal transition rates 
which are attained for large noise. 

In the non Markovian case the phase velocity, phase 
diffusion coefficient and Peclet number also prove phase 
synchronization between input and output. However the 
picture differs from the previous case showing a sequence 
of frequency locking modes. These different regions of 
locking are accompanied with low phase diffusion. The 
main conditions for locking are expressed by relations 
between the driving frequency fl and the time T spend 
on the excitation loop. The noise dependent and peri- 
odically modulated transition rates from the rest to the 
excited state act as a switch for the spiking. 1 : n lock- 



ing, i.e. a slow input and fast output, occurs for a cer- 
tain window of noise intensities if, in first approximation, 
the driving frequency is between 2ir/ ((n— 1/2)T) and 
27r/ (nT), respectively. For the opposite case of fast in- 
put and slow output we find 1 : n but also m : n, m < n 
frequency locking. The theoretical results for the non 
Markovian model of excitable systems agree well with 
simulations of a FitzHugh-Nagumo system. 

This work was supported by DFG-Sfb 555. The au- 
thors thank P. Talkner and J. A. Freund for fruitful col- 
laboration. 



Appendix A 

Our aim is to express the phase distribution V((f> — 
A(j),t — t) in terms of V(<l>, t) and its derivatives with 
respect to 4>, d n /d4> n V((f), t). To this end we start by 
expanding V(<p — A</>, t — r)) in a Taylor series around <fi 
and t, 

V{(t>- A&t-r) 

°° °° -i / an \ / am \ 

- £ E ass (s?) («= )»«.<><- A*n-rr 

n=0m=0 V V / \ / 

To process the time derivatives we use the Fokker-Planck 
equation eq. 1121 

taking care of the explicit time dependence of w(t) and 
D e g(t) which leads to 



v{<t>-A<f>,t-T)=v{<f>,t)- \a4>+ e 



(-r) m d m - x uj{t) 



I A( A 2 + v tit dm ' lD ^(t) 



dt T > 



Qfrn-l 

(-r) m d m ~ 2 uj(t) 
m(m-2)\ dt m ~ 2 



'WE 



-0(3). 



rn — 1 



(-r) m d" 1 " 1 ^) 
i! dt™- 1 



d 2 



where 0(3) denotes third or higher derivatives o{V{4>,t) 
with respect to 4>. The sums containing the derivatives 
of iv(t) and -D e ff(i) can be further evaluated, leading to 



and analogously 



y, (-T) m d^Desjt) 



m—1 



/ dr'D eff (t - t 1 ) 
Jo 



E 

m—1 



(-r) m d" 1 ' 1 ^) 
ml dt™- 1 



L^i ml dt m 1 



m—O 



and 



dr' u)(t — t ) 



E 



(-r) m d m - 2 uj{t) 
to(to-2)! dt m - 2 



dr'T'uj(t - r'). 
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Thus we eventually arrive at 
V(<f>- A<f>,t-r) = 



c) 



V{<)>,t) + (^'(t-r) - A4>)—V(ct>,t) + 



Next we insert our Ansatz eq. Ill -it 



oo fin 

p k (t) = ^qW(t) g-^t) 

n=0 V 



f>=27rfc 



-(A^ - A^ (t - r) + cT (i - r) — P(, , / . 
-0(3) 



(2), 



d 2 



where 



drD eff (r) + / dr(t - r)w(r). 
if 



into the dynamical equations eqs.© and JjJ 



P(i,*) = Ji^tPfc-i]- Ji^ 2 [Pfc] 

P(2,fc) = Jt^iPk] ~ Jt^iPk]- 



(Al) 
(A2) 



Using again the Fokker-Planck equation eq. ltl2)l the left 
hand side of eqs. i|Alj) and \A2l is given by 



A cx " 1 A f) 1 ^ f)n-\-l f) 



9^" 



n=0 1 n=0 r r 

The different terms on the right hand side of eqs. i|Alj) and l!A2j) read 

» f)n+l pin+2 fln+1 «n+2 

^[p,-!] - E^r^K^r - 2.— + ^_)p^ t) + ^ [c « q W](^- TT - 2*^)7>(<m) 



a.. , - 

+Jt j [cf ) q ( " ) ]^— t) + 0(n + 3) 



n+2 



n+2 



fin QTi-4- 1 /^Tl-(-2 



Equating now the coefficients of V(4>, t), -§gP{(j), t) an d 
TsV(^,t) finally leads to eqs. fl5 |t to l(T7 jl 
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